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\ The scalar contributions to the radiative decays of light vector 

\^ ' mesons into a pair of neutral pseudoscalars, V — > P^P^^, are studied 

■ within the framework of the Linear Sigma Model. This model has 
the advantage of incorporating not only the scalar resonances in an 

^ explicit way but also the constraints required by chiral symmetry. The 

Q_i' experimental data on (/> — > 7r'^7r'^7, (j) — > n^r]^, p -K^ir^'y and to 

^ , n^TT^'y are satisfactorily accommodated in our framework. Theoretical 

I predictions for (j) K^K^^, p T^^illi ^ ~^ ir^rj'y and the ratio 
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1 Introduction 



The radiative decays of the hght vector mesons {V = p,uj,(/)) into a pair of 
neutral pseudoscalars [P = 7T^,K^,ri), V —>■ P^P^'j, are an excellent labo- 
ratory for investigating the nature and extracting the properties of the light 
scalar meson resonances (5* = a, oq, fo). The reason is the following: looking 
at the quantum numbers of the initial vector and those of the final photon, 
both with J^^ = 1 , the system made of the two neutral pseudoscalars are 
mainly in a O^'^ state, i.e. a scalar state, or 2"*"+. However, the lightest tensor 
resonances have masses of the order of 1.2 GeV and therefore their contribu- 
tions to these processes are supposed to be negligible. In addition, there is 
also a vector meson contribution when one of the neutral pseudoscalars and 
the photon are produced by the exchange of an intermediate vector meson 
through the decay chain V — > VP^ —>■ P^P^'j. Fortunately, for most of the 
processes of interest the main contribution is by far the scalar one, thus mak- 
ing of the study of these radiative decays a very challenging subject in order 
to improve our knowledge on the lightest scalar mesons. This study also 
complements other analysis based on central production, D and J/ ip decays, 
etc. [1]. Particularly interesting are the so called golden processes, namely 
— >■ 7r''7r°7, TT^ri'y and p ir^n'^'y, which, as we will see, can provide 
us with valuable information on the properties of the /o(980), ao(980) and 
cr(600) resonances, respectively. A further motivation for the present work 
is the renewed interest on these radiative decays from both the theoretical 
[21 [31 H] and experimental [5], E] sides. 

At present, there are two experimental facilities that provide measure- 
ments on the V — > P^P^'j decays. One is the VEPP-2M e~^e~ collider in 
Novosibirsk with two experimental groups SND and CMD-2, and the other 
is the DA$NE 0-factory in Frascati with the KLOE Collaboration. The 
Russian experiment operates at different center of mass energies having the 
advantage of not only measuring the processes ir^n^'f and n'^ri'j 
but also p 7r°7r°7 and other similar decays. On the contrary, the Ital- 
ian experiment operates at a fixed center of mass energy around the mass 
and is only able, at least in principle, to measure the 0-decay processes but 
not others. However, due to its higher luminosity the statistical accuracy 
of DA$NE measurements is better than in VEPP-2M. Moreover, the good 
performance of the KLOE detector makes also the systematic error smaller 
and as a consequence the DA$NE measurements on 0-decay processes are in 
general more precise. For —>■ 7r°7r°7, the first measurements of this decay 
have been reported by the SND and CMD-2 Collaborations. For the branch- 
ing ratio they obtain 5(0 7r°7r°7) = (1.221 ± 0.098 ± 0.061) x 10"^ [7| 
and (0.92 ± 0.08 ± 0.06) x 10"^ [8], for m^^ > 700 MeV in the latter case. 
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More recently, the KLOE Collaboration has measured B{(f) n^TT^'-f) = 
(1.09 ± 0.03 ± 0.05) X 10"^ [9] in agreement with VEPP-2M results but 
with a considerably smaller error. In all the cases, the spectrum is clearly 
peaked at m,r7r — 970 MeV, as expected from an important /o(980) con- 
tribution. For 7r°?77, the branching ratios measured by the SND and 
CMD-2 Collaborations are 5(0 ^ 7r°r/7) = (8.8 ± 1.4 ± 0.9) x 10"^ pTO] 
and (9.0 ± 2.4 ± 1.0) x 10"^ [8J, and by the KLOE Collaboration are 5(0 ^ 
7r0r/7) = (8.51±0.51±0.57) x 10^^ from r/ 77 and (7.96±0.60±0.40) x 10"^ 
from rj —>■ 7r"'"7r~7r° [TT]. The two values are in agreement and also agree with 
those of VEPP-2M. Again, in all the cases, the observed invariant mass dis- 
tribution shows a significant enhancement at large n^r] invariant mass that 
is interpreted as a manifestation of the dominant contribution of the ao7 
intermediate state. For p —>■ 7r°7r°7, the only existing measurements in the 
literature come from the study of the e'^e~ — > 7r°7r'^7 process by the SND 
and CMD-2 experiments in the energy region 0.60-0.97 GeV. From the anal- 
ysis of the energy dependence of the measured cross section they obtain 
for the branching ratio B{p 7r°7r°7) = (4.llJ;^ ± 0.3) x 10"^ [12] and 
B{p 7r°7r°7) = (5.2+} J ± 0.6) x 10"^ [13], in agreement with the older 
result B{p ^ 7r°7r°7) = {A.StH ±0.2) x 10"^ [H]. These values can be 
explained by means of a significant contribution of the (77 intermediate state 
together with the well-known ujtt contribution. For uj —>■ 7r°7r°7, the val- 
ues of the branching ratio B{uj vr^vr 0^) = {6.6l:\i ± 0.6) X 10-5 [12] and 
B{uj 7r°7r°7) = (6.4+|^ ± 0.8) x 10"^ [lllj are found in agreement with the 
older GAMS result B{uj '^ 7r°7r°7) = (7.2±2.6) x 10"^ [15]. Finally, an upper 
limit for u — >■ 11^1]'-/ has been obtained: B{uj n'^ri'-f) < 3.3 x 10"^ at 90% 
CL [13j. For the rest of the processes, namely K^K^'-j and p ir^ri'y, 
the branching ratios are predicted to be very small and have not been yet 
measured. 

An early attempt to explain the V — > P^P^j decays was done in Ref. [16| 
using the vector meson dominance (VMD) model. In this framework, the 
V — *■ P^P^j decays proceed through the decay chain V —>■ VP^ -—>■ P^P^j, 
where the intermediate vector mesons exchanged are V = p for (cu, 0) —>■ 
ttOtt^, V = lu for p^ TT^n^-f, V = K*^,K*^ for -> K^K^-f, V = p,uj for 
(p, uj) — * 7r°?77, and V = p,u! for tt^tj'j. The calculated branching ratios 
are B^f^^'o = 1-2 x 10"^, BVMd ^ 5 4 ^ ^vmd = 1.1 x 10"^, 

= 2.8 X 10-^ 5y^°o^o = 2.7 X 10-12, = 4 X 10-10 and 

B^^^o^^ = 1-6 X 10-'' [16J. The first four are found to be substantially 
smaller than the experimental results quoted before. Later on, the same 
authors studied the V p^p^'y decays in a Chiral Perturbation Theory 
(ChPT) context enlarged to include on-shell vector mesons [IT]. In this for- 
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malism, 5J oo =5.1xlO-^5^ o =3.0x10-^ oo =9.5xl0-^ 
5^ = 2.1 X 10-^ ^0^0 = 7.6 X 10-^ 5^ = 3.9 X IQ-" and 
-B'^ =1.5x10^^. TakiiiEf into account these chiral contributions together 
with the former VMD contributions, one finally obtains bJ^^i^q = 6.1 x 

^Q_5 ^VMD+x ^ 3 g ^ ^VMD+x ^ 3.6 X 10-^ E™""?" ~ 5™ , 

oVMD+x _ r^x oVMD+x ^ T^VMD j oVMD+x _ t^VMD |^ 

for the first four, are still below the experimental results. Additional con- 
tributions are thus certainly required and the most natural candidates are 
the contributions coming from the exchange of scalar resonances (as stated 
in Sect. [H the contributions from tensor and higher spin resonances are neg- 
ligible). Needless to say that the previous two approaches do not contain 
the effect of scalar resonances in an explicit way. A first model including the 
scalar resonances explicitly is the no structure model, where the V —>■ P'^P^'j 
decays proceed through the decay chain V —>■ S'y —>■ P^P^^y, with S a scalar 
state, and the coupling VSj is considered as pointlike. This model seems to 
be rule out by experimental data on 7r°7r°7 decays [7j. A second model is 
the well-known kaon loop model (see Ref. [H] and references therein), where 
the initial vector decays into a pair of charged pseudoscalar mesons that af- 
ter the emission of a photon rescatter into a pair of neutral pseudoscalars 
through the exchange of scalar resonance^. In Ref. (TH], it was shown for 
the first time the convenience of studying the p^p^'j processes in order 
to investigate the nature of the /o and ao scalar mesons. In this pioneer 
work, B{(j) /o7 rnr-f) = 2.5 x 10'\ B{<p ao7 7r°r/7) = 2.0 x 10"^ 
and B{(j) —> (/q + 09)7 — * K'^K'^'y) = 1.3 x 10~®, for a four-quark structure 
of the scalar mesons involved, and 5.4 x 10~^, 2.4 x 10~^ and 2.0 x 10~^, 
respectively, for a two-quark structure. In addition, the contribution of the 
background processes is also calculated, B{(f) pvr ^ ini'-f) = 3.0 x 10^^ 
and B{(f) p°7r° — » T^^rjy) = 0.8 x 10"^. Due to its success in predict- 
ing the branching ratios as well as the mass spectra, the kaon loop model 
is used at VEPP-2M and DA$NE to extract information on the proper- 
ties and couplings of the /q and oq from the analysis of experimental data 
on ^ 7r°7r'^7 and (f) —>■ ir^ri'-f, respectively [7j-[ll]. A third model that 
makes use of the characteristics of the two former formalisms is developed 
in Refs. [20]-[26]. The previous three models include the scalar resonances 
ad hoc, i.e. the pseudoscalar rescattering amplitudes, where the scalar reso- 
nances are exchanged, are not chiral invariant. This problem is solved in the 
next two models which are based not only on the kaon loop model but also on 



-'^It is named the kaon loop model because it was first applied to the (j) 7r"7r"7 and 
— > 7r°7y7 processes [19]. 
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chiral symmetry. The first one is the Unitarized Chiral Perturbation Theory 
(UChPT) where the scalar resonances do not appear exphcitly but are gen- 
erated dynamically by unitarizing the one-loop pseudoscalar amplitudes. In 
this approach, taking into account only the contribution from unitarized chi- 
ral loops, = 8 X 10-^ B^l^'^^i = 8.7 X 10-^ B^^'^J « = 1.5 x 10-^ 

B^'^^l = 4.3 X 10-^ 5"^7o^o = 5 X 10-8, = 5.4 x 10"ii and 

B^^Jrn = ^ ^^'^ obtained 12HI [23 . A later analysis including 

more refined production mechanisms gives B^^J^^o^ = (1.2 ±0.3) x 10"^ and 

B^^J)^ = (0.6 ± 0.2) X 10"^ [30]. The second model is the Linear Sigma 
Model ^LcrM), a well-defined U{3) x U{3) chiral model which incorporates ab 
initio both the nonet of pseudoscalar mesons together with its chiral part- 
ner, the scalar mesons nonet. In this context, the scalar contributions to the 
V — s> P'^P'^'y decays are conveniently parametrized in terms of LaM ampli- 
tudes compatible with the corresponding ChPT rescattering amplitudes for 
low dimeson invariant masses. The advantage of the LcxM is to incorporate 
explicitly the effect of scalar meson poles while keeping the correct behaviour 
at low invariant masses expected from ChPT. 

The purpose of the paper is threefold. Firstly, to compute the (p — > 
K^K^j decay where the scalar effects are known to be dominant. This pro- 
cess is interesting to study, on one side, because it allows for a direct mea- 
surement of the KK couplings to the /o and mesons thus avoiding a model 
dependent extraction and, on the other side, since it could pose a background 
problem for testing CP-violation at DA<I>NE. The direct measurement of the 
couplings seems to be feasible in the near future with the higher luminosity 
expected at DAFNE-2. Having 50 fb~^, the number of expected K^K^'-f final 
state is in the range 2 ^ 8 x 10^ [5]. The analysis of CP-violating decays in 
— > K^K^ has been proposed as a way to measure the ratio e'/e [HI], but 
because this means looking for a very small effect, a 5(0 K^K^'-f) > 10~^ 
will limit the precision of such a measurement. Related to the — > K^K'^'j 
decay, there are also the processes — ^ (/o, 00)7 which are the main contri- 
butions to the former through the decay chain ^ (/o + cto)7 K^K^j. 
An accurate measurement of the production branching ratio and of the mass 
spectra for /o(980)/ao(980)7 decays can clarify the controversial nature 
of these well established scalar mesons. Secondly, to update our previous 
works on the decays — >• n'^ri'y ^32] and (p, a;) 7r°7r°7 [33] in view of the 
new experimental data provided by the KLOE [II], SND [12] and CMD-2 
|13j Collaborations. Concerning 7r°?77, the update consists in modifying 
the needed scalar amplitude to fulfill the chiral constraints at low energy, 
to incorporate the complete one-loop propagator, and to show the ex- 
plicit dependence of that new amplitude on the pseudoscalar mixing angle. 
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The intermediate vector meson exchange contribution is also included for 
comparison with the more precise forthcoming data. The new analysis of 
(p, uj) — s> 7r°7r°7 includes the most recent values of the mass and width of the 
a meson together with a determination of the size of the /o effects not consid- 
ered previously. Finally, to calculate for the fist time the scalar contributions 
to the p — i> 7r°r77 and u n^rj'y decays in the framework of the LcrM where 
the scalar meson effects are included in an explicit way. Our calculation will 
allow to quantify such contributions with the aim of establishing their rel- 
ative impact as compared to the vector meson exchange contributions. As 
a matter of completeness, we have also included in this work an updated 
version of our former analysis on the (j) — *• 7r°7r°7 decay [31] . 

The paper is divided as follows. In Sect. [21 we discuss and update the 
contributions from chiral loops in which our work is inspired. Sect. [3] is de- 
voted to the analysis of the needed scalar amplitudes in the framework of the 
LcrM, showing that these amplitudes can be considered as improved versions 
of their chiral counterparts. The VMD contributions that serve to complete 
our study are included in Sect. [H The final results where our predictions for 
the different V P^P^'-f decays are discussed in turn are given in Sect. [HI 
Finally, in Sect. [6] we present our conclusions and comments. A detailed cal- 
culation of the scalar amplitudes as well as the complete expressions of the 
one-loop scalar propagators and the invariant mass distributions are included 
in Appendix \M 



2 Chiral-loop predictions 

The vector meson initiated V — *• P^P^j decays cannot be treated in strict 
Chiral Perturbation Theory (ChPT). This theory has to be extended to in- 
corporate on-shell vector meson fields. At lowest order, this may be easily 
achieved by means of the 0{p'^) ChPT Lagrangian 

C2 = ^{D^U^D^U + M{U + U^)) , (1) 

where U = exp(i\/2P//), P is the usual pseudoscalar nonet matrixEl and, 
at this order, M = diag(m^, m^, 2m|- — ml) and f = = 92.4 MeV. The 
covariant derivative, now enlarged to include vector mesons, is defined as 
D^U = d^U - ieA^[Q, U] - ig[V^, U] with Q = diag(2/3, -1/3, -1/3) being 
the quark charge matrix and the additional matrix containing the nonet of 

^ In order to describe the processes involving the physical rj in the final state, the singlet 
term 770 has been added to the conventional octet part. 
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vector meson fields. We follow the conventional normalization for the vector 
nonet matrix such that the diagonal elements are (p^+ujq) / v^, (— p^+c^o) / a/2 
and 00) where luq = {uu + dd)/\/2 and 0o = ss stand for the ideally mixed 
states. The physical u and fields are approximately ujq — e0o and 0o + e^^o 
respectively, where e = ify — smipv — tan(^y = +0.059 ±0.004 accounts for 
the uj-(j) mixing angle in the flavour basis |35j . 

One easily observes that there is no tree-level contribution from the La- 
grangian ([T]) to the V P'^P^'-f amplitudes and that, at the one- loop level, 
one needs to compute the set of diagrams shown in Ref. [l7j. A straightfor- 
ward calculation of V{q*, e*) P'^{p)P^{p')'y{q, e) leads to finite amplitudes 
that are conveniently parametrized in the following way: 



, ^J, ._o^o 



-eg 



2\/27r 



2/m 2 



■m 



{a} L{mlo^o) X A 



X 

K+K-^ttOwO ' 



K+ 



(2) 



A, , 



= —A^ 



AX 



eg 



27r2m^+ 



{a}L(m^oi?o) X A\. 



(3) 
(4) 
(5) 



where {a} = (e* ■ e) {q* ■ q) — (e* ■ q) ■ q*) makes the amplitude Lorentz and 
gauge invariant, mpopo = s = {p + p'Y = {.<!* — is the invariant mass of 
the final dimeson system, and L{rn?popo) is the loop integral function defined 

as HZl [H [2HI EHl EZ! 



2{a~b) (a-6)2 
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and rj± 



1 ± Vl ^ 4z), a = my/rn?p+ and h = m?po / rn?p+ . The coupling 
constant g comes from the strong amphtude A{p tt+tt^) = —^/2ge* ■ 
iP+ ~ P~) with \g\ ~ 4.2 to agree with V{p vr+7r^)cxp = 150.3 MeV. 
However, for the decays we replace g by gg where \gs\ ~ 4.5 to agree with 
r(0 ^ ir+ir-)exp = 2.09 MeV [M] —in the good 5t/(3) hmit one should 
have \g\ = \gs\- These couplings are the part beyond standard ChPT which 
we have fixed phenomenologically. 

As seen from Eqs. (E]) and (HI), p ^ 7r°7r*^7 proceeds through a loop of 
charged kaons or pions, although the kaon contribution is suppressed by a 
factor of 10^ [T7]. (a;, 0) —>■ 7r°7r°7 proceed only by kaon loops since pion loop 
contributions are suppressed by G-parity for the lo case and, in addition, 
by the Zweig rule for the case — the same happens for the — >■ K^K^'y 
case, (p, 07,0) — s> 7r°?77 also proceed only by kaon loops because of isospin 
conservation. The four-pseudoscalar amplitudes are found to depend linearly 
on the variable s = ?7ipopo 

onlj{§: 
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C(pp 
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s — mt 



J TV 



(10) 



A^ 



'111 



wher^ fK = 1.22/^, 0p is the pseudoscalar mixing angle in the quark- 
flavour basi^ and (c0p,s0p) = (cos0p, sin0p). In this analysis we use 
0P = 41.8° {Op = —12.9°), a value obtained from the ratio — >• r]''y/r]'-y [38] 
and consistent with a fit to different decay channels |39j. It is worth noting 
that the amplitude in Eq. Qj includes not only the octet contribution (r/g) to 
the physical rj, as done in Ref. [17], but also the contribution from the singlet 
(fjo)- In doing so, we have enlarged the initial 5'f/(3)-flavour symmetry to 



TI16 ^p_^ p_ 



popo have not to be understood as four-pseudoscalar amplitudes in 

papa^ when computed using the 
chiral-loop framework of Ref. [17J. 

Strictly speaking, Jk — f-K at this order in the chiral expansion. 
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Table 1: Comparison between chiral- loop and LcrM predictions of branching 
ratios for different V P^P^'j decays. 

f/(3) — nonet symmetry — in order to obtain the relevant couplings for the 
singlet, thus incorporating rj-r]' mixing effects. A more general and rigorous 
extension of ChPT accounting for the effects of the pseudoscalar singlet in a 
perturbative way [40] would coincide at this order with our treatment of the 
mixing effects. 

Integrating the invariant mass distribution for the different V P^P^'j 
decays over the whole physical region one obtains the branching ratios shown 
in Table [H These results improve the predictions for these processes given 
in Ref. [U] where SU (3)-breaking effects were ignored and the pseudoscalar 
singlet contributions not studied. 

3 Scalar meson exchange 

We now turn to the contributions coming from scalar resonance exchange. 
From a ChPT perspective their effects are encoded in the low energy con- 
stants of the higher order pieces of the ChPT Lagrangian. However, the 
effects of the scalar states should manifest in the V p^p^'y decays not 
as a constant term but rather through more complex resonant amplitudes. 
In this section we propose scalar amplitudes which not only obey the ChPT 
dictates in the lowest part of the P^P^ spectra, as they must, but also gen- 
erate the scalar meson effects for the higher part of the spectra where the 
resonant poles should dominate. 

The Linear Sigma Model (LcxM) [ITl WI\ US] will be shown to be particu- 
larly appropriate for our purposes. In this context, the V P^P^'j decays 
proceed through a loop of charged pseudoscalar mesons emitted by the ini- 
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tial vector that, after the emission of a photon (in a gauge-invariant way), 
can rescatter into J^*^ = 0+"'" pairs of neutral pseudoscalars. The scalar 
contributions will be conveniently parametrized in terms of LaM amplitudes 
compatible with ChPT for low dimeson invariant masses. As an example of 
our procedure, we discuss in detail the K^K~ K^K^ rescattering ampli- 
tude needed for the K'^K^'j process. The remaining P^P~ popo 
rescattering amplitudes are written afterwards. 

The K^K^ K^K^ amplitude in the LaM turns out to be 



9R+R-R0R0 — 



gaR+R-daRORO 9 foR+ R' 9 foR° R° 



9aoR+R-9aoRORO 

s — m? 



^a^R^R0(0) 



s — m 



fo 



t 



'12) 



m 



ao 



where the various coupling constants are fixed within the model and can 
be expressed in terms of f^, /x, the masses of the pseudoscalar and scalar 
mesons involved in the process, and the scalar meson mixing angle in the 
flavour basis (f)s [HI US]- In particular, the 9k+r-roro coupling accounting 
for the constant four-pseudoscalar amplitude can be expressed in a more 
convenient form by imposing that the A{K^K~ K^)i^crM vanishes in the 

soft-pion(kaon) limit (either p — or p' — 0) — see Appendix |Al Then, the 
amplitude ( fT2l) can be rewritten as the sum of two terms each one depending 
only on s and t: 



A 



Lo-M 

R+R-^RORO 



m 



R 



^PR 



X 



""-^(c^.-v^s^^)^ 



R 



m 



D 



2P 



R 



(13) 



where Ds{s) are the S = a, /o,ao propagators — similar expressions holds 
for Daoit) — and (c05,s0s') = (costs', sin ^5). A Breit-Wigner propagator is 
used for the a, while for the /o and complete one- loop propagators are 
preferable (see App. E]) [T9|ll6]. 

A few remarks on the four-pseudoscalar amplitude in Eq. ( |T3i) and its 
comparison with the chiral-loop amplitude in Eq. ( fTTl) are of interest: 

i) for ms 



00 (S* = cr, /o, ao), the LcrM amplitude ( fT3|) reduces to 
s + t 



2ml 
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which agrees with the corresponding ChPT amphtude {^m^ — u)/2f'^ 
if the on-shell condition s + t + u = Am\ is invoked. This corresponds 
to the aforementioned complementarity between ChPT and the LcrM, 
thus making the whole analysis quite reliable. 

a) the L(tM amplitude ffT^ contains two terms: ^£o-m ^"^^ -^Lo-m- The 
former, being only s-dependent, can be directly introduced in Eq. ([5]) 
instead of ''^+k-^koko- Fortunately, this term contains most of the 
relevant dynamics due to the scalar resonances — mainly the /o(980) 
and ao(980) poles — for the process K^K^'^ where the dikaon mass 
spectrum covers the range Am\ < < "^i^. The contribution of the 

remaining term A^i^^yi^ being t-dependent, would require the calculation 
of a more involved four-propagator loop integral. However, this non- 
resonant ao contribution can be reasonably estimated by subtracting 
from the corresponding chiral-loop amplitude ffTTl) the resonant contri- 
butions from the a, /o and Oq terms in Eq. (fT^ taking mo-,/o,ao ~^ oo- 
The resulting expression corresponds to the desired non-resonant ao 
contribution in the iriaQ oo limit. The difference between this ao 
contribution and the same contribution for an ao of finite mass can be 
considered as negligible due to the high mass of the ao and the lack of 
poles in the t-channel. Therefore, an improved expression for the chiral- 
loop amplitude in Eq. (ITTl) taking into account the (pole dominated) 
s-channel scalar dynamics is the following: 



A 



Lo-M _ _ ^ 



4/. 



mi- — mi , I- -.0 



DJs) 



[C(ps - V/s 



R 

As desired, it has no t dependence and can thus be plugged in Eq. ([5]) 
instead of A^+^_^j^oko- Oi^s then obtains the s-dependent amplitude 



^^'^OKO-y — r,o^2 ^ ^^+R-'~^RORO , (15) 



which will be used from now on. Notice that for large scalar masses 
one recovers Eq. ([5]), that the /o(980), ao(980) (and, eventually, cr(600)) 
s-channel poles do now appear, and that the remaining term, (m|- — 
s/2)/2/|-, accounts for ao exchange effects in the t-channel. 
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Hi) the large widths of the scalar resonances break chiral symmetry if 
they are naively introduced in amplitudes, an effect already noticed 
in Ref. [17] . Accordingly, we introduce the /o(980) and ao(980) widths 
in the propagators only after chiral cancellation of constant terms in 
the amplitude (fT3l) . In this way the pseudo-Goldstone nature of kaons 
is preserved. 

iv) the K^K^ invariant mass spectrum for the K^K^'-y decay covers 
the region where the presence of the /o(980) and the ao(980) mesons 
should manifest. Because of the corresponding propagators in Eq. fll4l) . 
one should also be able to reproduce the effects of the /o and the ao 
poles at all K^K^ invariant mass values. 

Following the same procedure for the rest of relevant four-pseudoscalar 
amplitudes, one gets 



s — mt 



K 



mt^ - mt, 



m 



K 



A 



Lo-M 

K+K-^-K^ri 



s-m^m]^- mj^ 



2 



s/2 



(16) 



K 



m^ + m^ 



4/./. 



C(/)p 



C(f)p + 



A. 



Lo-M 

7r + 7r~ — ^TV^TV^ 



K 



2/2 

mt I mt 



(17) 



P 



m 



mt — m 



DJs) 



D 



As stated before, these amplitudes not only are compatible with their chiral- 
loop counterparts for low dimeson invariant masses but also generate the 
scalar poles that are dominant for higher invariant masses. The amplitudes 
f[T^ and (|T6l) - ([T8l) are a clear manifestation of the pursued complementarity 
between ChPT and the LcrM. However, Eq. ffT7|) is an exception. When all the 
scalar masses are sent to infinity the resulting amplitude does not match the 
chiral-loop one. This is because of the different treatment of S'[/(3)-breaking 
corrections in ChPT and the LaM. In Eq.(l9]), the Gell-Mann-Okubo mass 
relation Sm^^— 4m|-+m^ = has been used to express the amplitude in terms 
of m\ and m^, while in the LaM this relation is violated by second order 
S'C/(3)-breaking corrections and therefore only satisfied in the exact S'[/(3) 
limit. Consequently, the amplitudes ([9]) and ( JT71) can only be compared in the 
infinite scalar mass limit once the SU{'i) limit has been taken previously. In 
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that case, the pseudoscalar singlet is decoupled from the octet since Bp = 
(which implies (pp = arctanv^) and the whole pseudoscalar octet have a 
common mass, thus making both amplitudes to coincide in the ms oo 
limit. 

The final amplitudes Ay^^popo^ are obtained replacing ^p+p-^popo from 
the chiral-loop expressions in Eqs. ([2HS]) by the corresponding LcrM counter- 
parts A^p%._ ^popo shown above. Integrating the invariant mass spectrum 
for the different V P^P^'-f processes one finally obtains the branching ra- 
tios shown in Table [H As seen, taking into account the effect of scalar poles 
makes the LaM predictions bigger than the chiral-loop predictions. 



4 Vector exchange 

In addition to the LcrM contributions, which can be viewed as an improved 
version of the chiral-loop predictions, the analysis should be extended to 
include vector meson exchange in the t- and w-channel. These VMD contri- 
butions were already considered in Ref. p!6j. In this framework, — > P^P^'-f 
proceed through the exchange of intermediate (virtual) vector mesons mesons 
(y and V) in the direct and crossed channels of the V'^ —>■ {V, V^')P° —>■ 
P^P^'j decay chain. 

In order to describe these vector meson contributions we use the SU{3)- 
symmetric Lagrangians 

J^YYP = ^e^^'^'^id^Kd^V^P) , 

, (19) 

where G = is the upir coupling constant, / = and \g\ ~ 4.0 as 
follows from various p and uj decay data [35l HE] . The VMD amplitude for 
V{q*, e*) P°{p)P°{p')-f{q, e) is then found to be 

aVMD _ r(V^P^P^-AG^<' f PH'^}+{biP)} I P'\a}+{b{P')} \ (nn\ 
•Ay^popo^ - O^l/ 1- 1- ^M^-P^-iMvrv + M-^,-P2-iMv,rvJ ' ^^^^ 

with {a} as in Eq. ([2]) and a new amplitude 

{biP)} = -(e*.e)iq*-P)iq-P)-ie*.P)ie.P)iq*.q) 
+ (e* ■ g) (e ■ P) (g* • P) + (e • g*) (e* • P) (g ■ P) , 

where P = p + q and P' = p' + q are the momenta of the intermediate V 
and V mesons in the t- and u-channel, respectively. The intermediate vector 
mesons, V and V, can be either the u or the p mesons, with V = V in 
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7r°7r'^7 s-^id V ^ V in Tc^rj'-f] for (p K^K^'-f one obviously has V = K*^ 
and V = K*^ . The coefficient C is the same for both terms (using SU{'i)- 
symmetric couphngs) and changes from process to process according to well- 
known quark-model or nonet-symmetry rules: 

1 = C(p ^ 7r°7r°7) = 3C(a; ^ ttOtt^) = -^C{<P K'^K^-i) , (22) 

V2 

COS0P = 3C(p ^ vrV) = C'l^ ^ A7) , (23) 

and 

e = 3C(0 ^ 7r°7r°7) = sec 0p C(0 ^ tt V) , (24) 

for the (f) decays where the Zweig rule is operative. Notice the explicit de- 
pendence on the mixing angle for the processes involving the rj meson in the 
final state. Further details on these contributions are given in the following 
section. 



5 Results 

The final amplitudes A{y — P^P^'y) are thus the sum of the VMD contri- 
bution in Eq. fl20p plus the LcrM contribution containing the scalar resonance 
effects in Eq. (1141) and Eqs. (|TB1 - [TB|) . The P^P^ invariant mass distribution 
for these processes is given in Appendix IA.3[ The results of our analysis are 
now discuss in turn for each process. 

5.1 ^ K^K^-f 

As stated in the Introduction, this process is the only radiative decay where 
the /o and oq scalar mesons contribute simultaneously. This will allow, once 
this decay is measured presumably at DAFNE-2, to extract relevant infor- 
mation on the nature and couplings of both mesons, and to compare it with 
the one obtained from the already experimentally measured — n^Ti^'y and 
— i> n^ri'y decays. Therefore, a theoretical prediction for the different contri- 
butions to the branching ratio and mass spectrum of this process is welcome 
and useful. In addition, our analysis will serve to certify that the decay 
(f) — > K^K^'-y cannot pose a background problem for testing CP-violation at 
DA$NE since the calculated branching ratio is well below 10^^ (see below). 

The separate contributions from /o, and their interference, as well 
as the total result are shown in Fig. [H These mass spectra are computed 
assuming masses for the scalar resonances of m/g = 985 MeV and mag = 984.7 
MeV, and a pseudoscalar (scalar) mixing angle of 0p = 41.8° (0s = —8°). 
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Figure 1: dB{(j) ^ K'^ K'^^)/ dm j^oko x 10^ (in units of MeV~^) as a function of 
the dikaon invariant mass rrij^oj^o (in MeV). The dotted, dashed and dot-dashed 
hues correspond to the separate contributions from /o, oq and their interference, 
respectively. The solid line is the total result. 



The values of the /o mass and the scalar mixing angle are obtained from the 
analysis (see the dedicated subsection) whereas the oq mass is 
taken from Ref. [35] and (pp from the ratio — > i]''~f/ri'~f [38j. As seen, the 
/o contributes more strongly than the due to a smaller imaginary part of 
the propagator and a larger coupling to kaons. The interference is negative 
since isospin invariance implies Qf^K+K- = QfoK^R'^ and gaoK+K- = ~9aoK°K°- 
Integrating the K^K^ invariant mass spectrum one obtains for the scalar 
contribution B{(j) — ^ K^K^'-))i^(jm = 7.5 x 10~^. It is worth mentioning that 
the value obtained is very sensitive to the scalar mixing angle, a change of 
1° modifies the branching ratio around 20%. The whole scalar contribution 
includes not only the /o and oq effects but also the a ones. However, the 
latter are negligible due to the suppression of the aKK coupling i^ mg. ~ 
and for kinematical reasons. Numerically, they amount to less than 1% of 
the total result. Finally, the exchange of intermediate K*^ and K*^ vector 
mesons is calculated to be 5(0 K°K^-f)yMD = 2.0 x lO'^^ ^nd therefore 
negligible. In summary, our predicted branching ratio including the scalar 
resonances explicitly is 5(0 — > K^K^'j) = 7.5 x 10~®. This value is in 
agreement with previous phenomenological estimates [271 ESI EO] (see also 

^ A value of rricr — 478 MeV is taken from the analysis of £)+ 7r^7r+7r+ performed 
by the E791 Coll. [49] at Fermilab. 
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Ref. [Sn] for a review of earlier predictions). Notice that the branching ratio 
obtained here is one order of magnitude larger than the chiral-loop prediction 
B{(j) — s> K^K^'-j) = 4.1 X 10~^. However, it is still one order of magnitude 
smaller than the limit, O{10~^), in order to pose a background problem for 
testing CP-violating decays at DA$NE. 



5.2 (J)^ /o7/ao7 

In the kaon loop model, these two processes are driven by the decay chain 
(f) —>■ K~^K~{'y) /o7 and ao7. The amplitudes are given by 

^ = 2n^m^ ^""^ ^("^/o(ao)) X 9fo(ao)K+K- , (25) 
where the scalar coupling constants are fixed within the LaM to 

9foK+K- = TTi y^^S + V2C05) , gaoK+K- = TT? ■ (26) 

^JK ^JK 

The ratio of the two branching ratios is thus 

- (1 - ml/,nir " - + • <2^) 

where the approximation is valid for mj^ maQ. For (ps = —8°, one gets 
R^^Z^gy/ao-y — 1-6 which should be compared with the experimental value 
R^l^fy/ao-y = 6.1 ± 0.6 [llj. However, this value is obtained from a large 
destructive interference between the /o7 and a'j contributions to ^ n^n^'j, 
in disagreement with other experiments [7]. In addition, the approximate 
expression for the ratio (1271) is only valid when the /o mass is below the 
charged kaon threshold (2m^+ ~ 987 MeV). If not, the steep behaviour of 
the loop function after threshold makes the approximation meaningless and 
the exact expression has to be taken. In this case, B^'^^^/^^^ can be much 
smaller than the given prediction (independently of the mixing angle value). 
Conversely, if the mass, which we have kept fixed to rriao ~ 984.7 MeV, 
is moved to a value above threshold, then the ratio can be even larger than 
the experimental result. Therefore, a precise fit of the /o and masses 
using the n^ir^'y and n^rj'y decays is mandatory before drawing definite 
conclusions on R^^f^^/aQ-y- Once these masses are fixed and assuming the 
validity of the kaon loop model as the production mechanism and the LaM 
as the model for fixing the SPP couplings, the measurement of this ratio 
could be used to get some insight into the value of the scalar mixing angle. 
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Figure 2: dB{(f) 7r°7r°7)/dm^o^o x 10® (in units of MeV~^) as a function of the 
dipion invariant mass m^o^o (in MeV). The dotted, dashed and dot-dashed hues 
correspond to the separate contributions from VMD, LuM and their interference, 
respectively. The sohd line is the total result. The long-dashed line is the chiral- 
loop prediction. Experimental data are taken from Ref. [7J (solid star) and Ref. [9] 
(open diamond). 

5.3 (j) 7rV7 

This process is considered the first of the so-called golden processes since the 
/o contribution is seen to be dominant. Hence, the precise measurements of 
the total branching ratio and mass spectrum performed at SND [7], CMD-2 
[8] and KLOE [9] allow to extract valuable information on the properties and 
couplings of this resonance. In this analysis, we predict the scalar effects on 
this interesting and important process using the LaM as the framework for 
calculating the /o and a resonant contributions. In addition, we also predict 
the vector meson exchange effects. 

The 7r°7r° invariant mass distribution, with the separate contributions 
from the LcrM, VMD and their interference, as well as the total result, are 
shown in Fig. El We use = 478 MeV [19], = 258 MeV, as required by 
the LcrM, = 985 MeV and 05 = —8°. The latter two values are obtained 
from the best fit to experimental data in Refs. [TJ E] and are the only two 
free parameters of our prediction. The global agreement with the data is 
rather acceptable. As expected, /o(980) scalar meson exchange contributes 
decisively to achieve this agreement. Indeed, for the integrated branching 
ratio one obtains 5(0 vr°7r°7)Lo-M+VMD = l-H x 10""^, quite in line with the 
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experimental results quoted in Refs. [3, E] and previous analyses [IHl [211 ED] • 
Without the VMD contribution, our result decreases only by some 10%, 
B{(f) — » vr'^7r''7)Lo-M = 1-00 x 10"^, but it still remains well above the chiral- 
loop value, B{(f) ti^tt^^)^ = 4.2 x 10~^, which did not contain the crucial 
scalar contributions. 

The reason for not achieving a better accord between our prediction for 
the TT^vr'^ mass spectrum and the KLOE experimental data is the simplicity 
of our signal amplitude. It is derived from a tree level LaM where the scalar 
contributions to the s-channel, the dominant production mechanism, have 
been kept explicitly while other contributions, scalar or not, have been inte- 
grated out in a manner that preserves the chiral-loop result (see App. lA.ll for 
details). The advantage of using such a simple amplitude is its predictivity, 
i.e. with only two parameters (the couplings of a and /o to pions and kaons 
are fixed within the model and the a mass is fixed from experiment) one is 
able to get a reasonable description of the mass spectrum and a better value 
for the integrated branching ratio. Besides, our background amplitude is de- 
duced from a VMD model where the required 0p7r° and p7r°7 couplings are 
written in terms of a single parameter g which is also fixed from experiment. 
A way of improving our amplitude could be on the side of the background to 
include an additional phase to take into account pn scattering effects and fix 
(ppn^ and pTc^j independently from experiment. On the side of the signal, it 
would be welcome to incorporate next-to-leading effects in the mixing ma- 
trix of /o and a, beyond the tree level description in terms of a mixing angle, 
and include more decay channels in the propagators. In this direction, an 
interesting analysis is that of Ref. [1] where a more involved amplitude also 
inspired by the LcxM is used for describing simultaneously 7r°7r'^7 and 
TTTT scattering. Their amplitudes accommodate all the above improvements 
as well as a phase for the elastic background of the vrvr and KK scattering. 
As a result, the agreement between their predictions and the KLOE experi- 
mental data is significantly better than ours. Nevertheless, the values for the 
scalar masses and the a width they obtained from their fits are close to ours 
(see Ref. |1]) and the hierarchy of the scalar couplings follows the pattern 
dictated by the LaM, thus giving more consistency to our simpler approach. 

It is worth noting that our prediction is very sensitive to the /o mass and 
scalar mixing angle. As an example, B{(j) n^ir'^'y) = (1.20,0.98) x 10~^ 
for = (980,990) MeV and B{<p 7rV7) = (1.16,1.06) x 10"^ for 
(ps = (—9°, —7°). The dependence on both parameters is because of the /ovrvr 
coupling which in the LcrM is proportional to mj^ and sin05. This coupling 
enters into the calculation not only in the scalar part of the amplitude but also 
on the imaginary part of the propagator (remember that we use a complete 
one- loop propagator for this resonance). In addition, if the /o mass is above 
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Figure 3: dB{(j) 

^'^0^0 (ill GeV^) — horizontal axis — and w-^o^ (in GeV^) — vertical axis — . The 
left, centered and right Dahtz plots correspond to the VMD, LuM and total con- 
tributions, respectively. 



the two charged kaon threshold the kaon-loop function drops drastically thus 
decreasing the /o contribution. Consequently, the cf) 7r°7r°7 decay could 
be used to obtain relevant information on these two parameters. Concerning 
the a contribution, it is suppressed with respect to the /o, Smce (JcrKK 
(m^ — m\) ~ for m„ ~ rriK, but it still amounts to 10% of the total result. 
By contrast, the chiral-loop prediction shows no suppression in the region 
m,r7r — 500 MeV, see Fig. O The effects of the a are enhanced if instead of 
using rrifj = 478 MeV, one takes rricr = 513 MeV from the Dalitz analysis of 
Klir+Tf- by the CLEO Coll. In this case, however, the fit to the 
^0^0 j^g^gg spectrum is worse and the branching ratio incompatible with the 
KLOE result. On the contrary, our prediction is hardly sensitive to the a 
decay width. Changing from the LcrM value Fo- = 258 MeV for mg- = 478 
MeV to F^ = 324 MeV [4.9J or F^ = 335 MeV [5l] makes no substantial 
difference (less than 1% in both cases). Therefore, the dependence on the a 
parameters, particularly on m^, could be used in the near future to fit its 
value once more precise experimental data on — > 7r''7r''7 is available. Finally, 
the VMD contribution is calculated to be 5(0 vr*^7r°7)vMD = 8.3 x 10~^ 
that although small should not be considered as negligible since as stated 
before it is responsible together with its interference with the scalar part 
for a 10% of the final result. For the sake of comparison with experiment, 
we also show in Fig. [H] the predicted Dahtz plots for the scalar and vector 
meson exchange contributions, as well as the total result. As seen, the VMD 
plot shows the exchange of intermediate p mesons at tti^o^ — nT-f, (horizontal 
band) and tti^o^ + '^^o^o — 2m^o + fn"^ — fn"^ (inclined band), while the LaM 
plot shows the exchange of an /o(980) at ttl'^o^o ~ (vertical band). 
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One way to further test the goodness of our approach apphed to this 
decay, in particular for the amphtude used, would be to perform a simulta- 
neous analysis of — ^ and TTTT scattering along the lines of Ref. [52] . 
By unitarity, the phase associated to the / = J = partial wave of the 
K^K~ — >■ tt^tt'^ amplitude must be equal below the kaon threshold to the 
phase dQijiTx — » txtx). This kind of analysis, which is beyond the scope of 
the present one and left for future work, has been already carried out for 
the SND data in Ref. [18] and for the more precise KLOE data in Ref. 
showing in both cases that with an amplitude more elaborate than ours (see 
previous discussion) it is possible to get simultaneously a remarkable fit to 
the TT^TT^ mass spectrum and the vrvr phaseshift. 

5.4 ^ 7r^?77 

The second of the golden processes, cj) ir'^ri'j, is interesting to analyze since 
it provides relevant information on the properties and couplings of the Oq 
resonance, whose effects on the process are also seen to be experimentally 
dominant. Here, we present a parameter- free prediction for the different 
contributions, scalar and vector, to this decay using the LcxM for calculating 
the KK irrj invariant amplitude. Our calculation of the scalar contribution 
not only satisfies the chiral constraints at low energy but also incorporates 
the pole effects of the ag in an explicit way. In this manner, we are able to 
nicely reproduce the peak and the lower tail of the mass spectrum at large 
and small ttt] invariant masses, respectively. 

The separate contributions from the LcrM, VMD and their interference, 
as well as the total result, are displayed in Fig. HI We set rriaQ = 984.7 
MeV [33] and (pp = 41.8° [HH]. As seen, the agreement with data is pretty 
good and the Qq effects dominate. The integrated branching ratio is B{(f) 
7r°r77)LcrM+VMD = 8.2 x 10~^, in accord with the experimental values in 
Refs. [ini HI] and previous estimates [HI ESI EOl E3]. The sole contribu- 
tion of the ao(980) amounts to B{(f) —> 7r*^?77)LcrM = 7.8 x 10~^, a value three 
times bigger than the chiral-loop prediction B{(f) —>■ Tc^rj'-f)^ = 2.9 x 10~^. 
The contribution from vector exchange — mainly p exchange since the uj ef- 
fects are calculated to be less than 1% of the total value — is found to be 
B{(j) — i> 7r°?77)vMD = 3.4 x 10~^, a small effect, around 5% of the final result 
including its interference with the scalar part, which however should be taken 
into account in future analyses with more precise experimental data. The 
obtained branching ratio is rather sensitive to the T]-ri' mixing angle. If one 
uses 6p = arcsin(— 1/3) ~ —19.5° (0p ~ 35.3°) the new result differs in a 
10%, a sizable effect which is however of the same order of the current exper- 
imental errors. It is worth remarking that our predictions for the invariant 
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Figure 4: dB{(j) 7r^r]j)/dm^o^ x 10'' (in MeV^-*^) versus m^o^ (in MeV). The 
curves follow the same conventions as in Fig. [2j Experimental data are taken 
from Ref. [10] (solid star) and Ref. [11]: (open diamond) from ^ 77 and (solid 
diamond) from r] — > tt~^it~tt^. 

mass distribution and branching ratio of the (j) — > n^rj'j decay are free from 
adjustable parameters, i.e. the ao mass and the pseudoscalar mixing angle 
are fixed while the Oq couplings are predicted within the LaM. A main draw- 
back of this approach is the large value of the aoTTi] coupling constant which 
implies F^^^^^ = 386 MeV. However, this large value is considerably reduced 
once the well-known Flatte corrections for the Oq propagator are introduced 
|54j . The complete one-loop propagator used in this analysis is nothing else 
than an improved version of the Flatte formula taken into account the full 
energy dependence. 

In summary, our present analysis supersedes Ref. [32] incorporating not 
only the chiral constraints at low energy of the (p 7r°?77 amplitude but also 
the explicit dependence on the rj-rj' mixing angle and the complete ao prop- 
agator in the corresponding four-pseudoscalar amplitude, three new features 
which were not taken into account in the former work. 

5.5 p TT^TT^^ 

The last of the golden processes is the p -K^ir^'-f decay, where the o"(600) 
contribution is supposed to dominate the scalar part of the amplitude since 
the /o(980) effects are suppressed by kinematics. In addition, there is also 
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Figure 5: dB{p 7r°7r°7)/dm^o^o x 10'' (in MeV"^) versus m^o^o (in MeV). The 
dotted, dashed and dot-dashed hnes correspond to the separate contributions from 
VMD, L(tM and their interference, respectively. The sohd line is the total result. 

a well-known intermediate utt contribution which accounts for an important 
part of the total signal. In spite of this sizable vector contribution, the 
p — > 7r°7r°7 decay should be considered as an ideal process for extracting the 
mass and decay width of the controversial and not well established a meson. 
Therefore, a precise measurement of the branching ratio (already existing) 
together with its invariant mass distribution (not yet available) could help to 
fix unambiguously the a properties and discriminate among different models 
predicting the size and shape of its contribution. In this analysis, we provide 
a theoretical prediction for the branching ratio and mass spectrum of the 
p TT^TT^'-f decay taking into account the contributions of the a and /o 
scalar mesons, calculated within the framework of the LcrM, as well as the 
UJTT vector contribution given by the VMD model. 

The separate contributions from the LcxM, VMD and their interference, 
as well as the total result, are displayed in Fig. We use m„ = 478 MeV 
[19], = 258 MeV, as required by the LcxM, m/^ = 985 MeV and 0^ = -8°. 
The latter two values are fitted from the comparison to experimental data of 
our prediction for — > ir^n^'j (see Sect. 15.31 for details). We have chosen a 
simple Breit-Wigner propagator for the a resonance. This parametrization is 
considered in the literature as the preferred choice when dealing with broad 
scalar states. As seen from Fig. the scalar contribution is peaked at m^^r — 
500 MeV, as expected from a single resonance of that mass which is far from 
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thresholds, and the interference with the vector contribution is positive in the 
whole mass range and suppressed in that particular region. The integrated 
branching ratio obtained, B{p 7r°7r°7)Lo-M+VMD = 4.2 x 10^^, is in perfect 
agreement with the new experimental results [I2l [13] (confirming the older 
measurement [H]) and previous analyses [211 EE], EH [29l [55] . The scalar part 
alone amounts to B{p vr°7r°7)LcrM = 2.2 x 10~^, which is clearly dominated 
by the a contribution — /q effects are found to be less than 1% of the total 
result — , while the intermediate ujtt^ exchange leads to B{p 7r°7r°7)vMD = 
9.0 X 10-6. Notice from the comparison of the former scalar contribution 
with the chiral-loop prediction, B{p —>■ TT^Tc^'-f)^ = 1.1 x 10"^, that taking 
into account the a effects in an explicit way is mandatory in order to agree 
with experimental data. If not, a lower less preferred value is obtained. The 
cr contribution is slightly sensitive to the scalar mixing angle since (ps — ~8° 
and gamr oc COSTS' ~ 1 for small values of 0^. However, it is quite dependent 
on the a properties, mass and width. To study this dependence, we have 
redone our analysis applying different values for m„ and T^. With the same 
mass value, rria- = 478 MeV, but replacing the predicted LcrM width by 
= 324 MeV [IS] one gets B{p 7r°7r°7)LaM = 1-5 x 10-^ which implies 
B{p 7r°7r°7)Lo-M+VMD = 3.3 x 10"^. If instead, m^j = 513 MeV and 
To- = 335 MeV are taken from Ref. [HI], one gets B{p 7r°7r°7)Lo-M = 

1.5 X 10"^ and B{p vr°7r°7)L<TM+VMD = 3.5 x 10~^. Finally, the chiral- 
loop prediction, which is formally equivalent to the mo- — > oo limit, gives 
B{p 7r°7r°7)vMD+x ~ ^-^ ^ 10~^. The various predictions for the 7r°7r° 
mass spectrum using the former values of a mass and width are shown in 
Fig. [6l As seen, a detailed experimental check of this mass spectra together 
with a more precise measurement of the branching ratio would help to fix 
the mass and width of the elusive a meson. 

This new analysis, taking into account the latest values of the Breit- 
Wigner mass and width of the a resonance as well as the effects of the /o 
contribution for the first time, allow us to update our former predictions for 
the shape and branching ratio of the p —>■ n^n^'j decay [33j and confront 
them with the most recent experimental data. 

5.6 u 7r^7r^7 

This decay is less interesting from the point of view of extracting the prop- 
erties of the a meson. The reason is the following. Because of G-parity pion 
loops are forbidden and the chiral scalar contribution to the process is then 
driven by kaon loops. Therefore, due to the large kaon mass, it is expected 
to be small as compared to the p — > ir^n^'j case where pion loops are allowed. 
In consequence, the u —>■ n^Tc^'-f transition is dominated by the VMD con- 
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Figure 6: dB{p tt^tt °7)/dm^o^o X 10^ (MeV~^) versus m^o^o (MeV). The 
various predictions are for the input values: mo- = 478 MeV and To- = 324 MeV 
[19] (sohd hne); = 478 MeV [l9] and = 258 MeV, as predicted by the LcjM 
(dot-dashed line); and nio- = 513 MeV and To- = 335 MeV [51], (dashed line). The 
chiral-loop prediction is also included for comparison (dotted line). 

tribution. Recently, Guetta and Singer |56j have explored the possibility of 
incorporating the pion loops through p-u mixing effects. However, as we see 
below, these effects are not enough to make the scalar contribution relevant. 

Ignoring for the moment p-uj mixing, i.e. assuming that the physical 
u = u)^^^ with no / = 1 contaminations, the well-known p meson ex- 
change contribution via the VMD decay chain u) —>■ pir^ —>■ 7r°7r°7 gives 
T{uj 7r°7r°7)vMD = 244 eV and B{uj 7r°7r°7)vMD = 2.9 x 10-^ a value 
which is more than two standard deviations below the experimental mea- 
surement in Ref . [1^ . This result is improved using a momentum dependent 
width for the p meson. 



leading to T{u — > vr°7r°7)vMD = 272 eV, which is some 12% larger than 
the previous value, as already noticed in Ref. [56]. There is also another 
contribution to the u 7r*^7r^7 amplitude coming from chiral loops. However, 
as stated before, this chiral loop contribution (given only by kaon loops in 
the good isospin limit with uj = uj'^^) is very small, r(a; — ^ ir^TT^'-f)^ = 1.3 
eV, and can be safely neglected. Its improved version taking into account a 



rp(g^) 




(28) 
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resonance effects is also negligible, T{lj 7r°7r°7)Lo-M = 1-5 eV, due to the 
suppression of the g^rKK coupling for m„ ~ uik (see Sec. 15.31 for details). 

Taking now into account p-u mixing, there is an indirect contribution 
to the u TT^n^j amplitude that is written as A^^^{u! — vr°7r°7)vMD + 
eA{p 7r°7r°7), where A^=^{uj 7r°7r°7)vMD = ^ vr°7r°7)vMD —the 

proportionality factor follows from the SU (3)-symmetric Lagrangian (1191) — 
and e is the p-tu mixing parameter given bju 



-0.003 + i 0.034 , (30) 



with Ml^{ml) = (-3800 ±370) MeV^ [57J. Approximating the new, isospm 
violating term of u ^ 7r°7r°7 by the VMD contribution €A{p 7r°7r°7)vMD5 
one increases the previous estimate to r(co' ^ tt^tc^'j) = 296 eV. However, a 
more complete treatment, with A{uj — > n^n^-y) = A^^'^^uj —>■ 7r°7r*^7)vMD + 
eA{p vr°7r°7)vMD+Lo-M, seems preferable. The 7r°7r° invariant mass spectra 
corresponding to this amplitude have been calculated with the same rricr and 
Tcr values introduced for the p 7r°7r°7 case. The sensitivity on the cr 
parameters is seen to be minimal and all the results almost coincide with the 
curve for = 478 MeV and = 324 MeV [49j plotted in Fig. El 

The integrated width and branching ratio for the former values are pre- 
dicted to be r(co' ^ 7r°7r°7)vMD+LaM = 297 eV and B{iu ^ 7r°7r°7)vMD+LCTM = 

3.5 X 10~^. If all a meson effects are neglected — chiral-loop prediction — one 
then obtains r(a; —>■ 7r°7r°7)vMD+x = 307 eV and B{u — > vr°7r°7)vMD+x = 

3.6 X 10~^, only a 4% above the previous results and hardly distinguish- 
able (see Fig. [7]). Notice that our values are still substantially lower than the 
central value reported in Ref. [12]. The reason is that we are using an SU (3)- 
symmetric formalism where all the VVP and VP7 couplings are deduced from 
the single VPP coupling g taken from r(p 7t~^tt^). In this formalism, the 
product of couplings guipTrdp^TrO-y = = 3.68 GeV~^ for \g\ = 4.24 and 
g^p^ = G = 14.8 GeV^^. Alternatively, one could use the experimental val- 
ues g^^^ = 16.7 GeV^^ and fi'^o^o^ = 0.25 GeV~^, extracted from the study 
of the process e~^e~ — > ujti^ 7r^7r°7 above 1 GeV at CMD-2 [5H] and the 
current value of r(p° 7r°7) = (91 ± 13) keV [35], to fix gZ^f^g^p^% = 4.24 



^ A further effect of this p-uj mixing is to replace the p propagator in by 



1 + ^""^ ' I , (29) 



7W 



Dp{s) D<p{s) Dp{s) \ gp^^ D^{s) _ 

with Dv{s) = s — rriy + i myFy for V — p,uj and gu-K-y/gp-Ki = 3 in the SU (3)-symmetric 
VMD framework. 
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Figure 7: dB(uj tt^tt °'y)/dm^o^o X lO'^ (MeV"^) versus m^o^o (MeV). The 
plotted curves are for mg- = 478 MeV and To- = 324 MeV (solid line) and dropping 
all a meson effects — chiral loop prediction — (dotted line). 

GeV~^, which is 15% larger than the 5'f/(3)-symmetric result. This enhance- 
ment applied to the squared amplitude gives T{lj — > 7r°7r°7)vMD+Lo-M = 394 
eV and B{u — > vr°7r°7)vMD+Lo-M = 4.6 x 10~^, a value which is still below 
the experimental result in Ref. [12] by more than one standard deviation. 
Our prediction is in fair agreement with a previous analysis in Ref. [22] but 
clearly disagrees with Ref. [20]. 

In conclusion, our final results for the branching ratio of the a; 7r°7r°7 
decay is B{uj vr°7r°7)vMD+L(TM = (3.5-4.6) x 10~^, where the numeric in- 
terval takes into account the departure from the S'?7(3)-symmetric VMD 
prediction. Therefore, 4.6 x 10^^ seems to be the maximum branching ra- 
tio acceptable within VMD model and additional theoretical study would be 
required to explain the large value of B{uj n'^n^'-f). 

5.7 p 7T^r]j and u 7r^?77 

The scalar contribution to these two processes is driven through a loop of 
charged kaons since the J = 1 of the final state tit] forbids a loop of pions due 
to Bose symmetry. Therefore, scalar effects are expected to be suppressed, 
because of the large kaon mass, and vector meson exchange contributions 
should dominate. Consequently, these decays seem to be not very promising 
from the point of view of extracting the ao(980) properties. However, they 
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Figure 8: dB{p ^ 7r°7?7)/dm^o^ x 10^^ (MeV-^) —left plot— and dB{u; 
'K^rj'^) / dm^o^ X 10^ (MeV~-^) — right plot — versus m^o^ (MeV). The curves follow 
the same conventions as in Fig. \5\ 

can still be used to test the hypothesis of a dominant intermediate vector 
contribution as soon as experimental data is available. Unfortunately, the 
size of the predicted branching ratios, of the order of 10~^ for p TT^rj'y 
and 10~^ for lj n^ri'y, make of their experimental analyses a very difficult 
task not feasible in the near future. At the moment, only the upper limit 
B{uj 7r°r/7) < 3.3 x 10"^ at 90% CL has been obtained [13]. 

The contributions to p — 7r°?77 and u — > n^ri'j from the LcrM, VMD 
and their interference, as well as the total result, are shown on the left and 
right hand sides of Fig. [HI respectively. We have used the same input values 
for the oq mass and pseudoscalar mixing angle as in the — Tr^rj'y anal- 
ysis. As seen, the vector exchange contribution is dominant in both cases 
and the integrated branching ratio gives B{p —>■ 7r''?77)vMD = 4.5 x 10~^° 
and B{uj vr*^r77)vMD = 9.9 x 10~®. The difference of two orders of 
magnitude in branching ratio is easily understood in terms of the ratio 
^VMD /^VMD ^ T^/(9To) where the numerical factor is deduced from 
Eq. fl^jl) . Although smaller in size, the scalar contribution to the p — > 
7r°r77 decay is more relevant, as compared to the vector one, than for the 
uo TT^'q'y case. In the former, its integrated branching ratio amounts to 
B{p — > 7r°?77)Lo-M = 1-3 X 10~^°, around 25% of the total signal, while 
in the latter B{uj '^^f]l)LaM = 3.4 x 10~^, less than 5% of the signal. 
If the Oq resonance is removed from the analysis, the prediction of chiral 
loops is B{p Ti^rj'^)^ = 6.3 x lO"^"*^ and B{uj —>■ ri'^rj'^)^ = 1.6 x 10~^, 
so in both factor of two of enhancement is achieved when the oq 

is introduced explicitly. In view of all the above discussion, it would be 
preferable as long as experimentally possible to look for the Oq effects in 
the p — > n^rij decay since in the uj — > 7r''?77 they are not discernible (see 
Fig. M for comparison). However, as stated before, p n^ri'-f should be 
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Figure 9: dB{p ^ 7r°r/7)/(im^o^ x 10^^ (MeV-^) —left plot— and dB{oj 
Tr^r]^)/dm^o^ x 10^ (MeV~-^) — right plot — versus m^o^ (MeV). The plotted curves 
are for an explicit qq resonance with mag = 984.7 MeV (solid line) and dropping 
all uq meson effects — chiral loop prediction — (dashed line). 

considered as a very rare decay and hence very hard to detect. Fortunately, 
one should keep in mind that there is a golden process for extracting the oq 
properties, which is obviously the — >• 7r°?77 decay discussed previously. Fi- 
nally, our final results including scalar and vector exchange contributions are 
B{p 7r°r/7)vMD+LaM = 5.2x10"^° and5(^ T^°Vl)vMD+LaM = 9.7x10"^ 
which should be compared with B{p — > 7r°777)vMD+x = 4.7 x 10^^° and 
B{lj 7r''?77)vMD+x = 9-7 ^ 10"*^ in absence of scalar effects. 

6 Conclusions 

In this work we have performed an extensive and exhaustive analysis of the 
scalar and vector meson exchange contributions to the V —>■ P^P^'y decays 
with V = p,uj,(j) and P^P^ = h^tt^ ^Tc^rj, K^K^ . The scalar contributions, 
which turn out to be dominant in most of the interesting cases, have been 
studied within the framework of the LaM. This model has the advantage 
of incorporating not only the scalar resonances in an explicit way but also 
the constraints required by chiral symmetry. The complementarity between 
ChPT and the LaM have been used to parametrize the needed scalar am- 
plitudes making explicitly the effect of scalar meson poles while keeping the 
correct behaviour at low invariant masses. Therefore, our treatment of the 
scalar contributions can be considered as an improvement of the correspond- 
ing chiral-loop predictions. Moreover, we have updated the vector exchange 
contributions. A summary of the results obtained is presented in Table [21 
Our theoretical predictions are all compatible with experimental data. This 
nice agreement is achieved not only for the branching ratios but also for 
the available mass spectra. In general, we have shown that (p tt^tt^'j, 
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S(y^pOp%) LaU VMD L(7M+VMD 

^ 7rV7 1.00 X 10-^ 8.3 x 10-^ 1.11 x 10"^ 
(j) 7r°r/7 7.8 x 10"^ 3.4 x 10"^ 8.2 x 10"^ 

^ K^K^-f 7.5 X 10-^ 2.0 X 10"^^ 7 5 ^ iQ-s 

p ^ 7r°7r°7 2.2 x 10"^ 9.0 x 10"^ 4.2 x 10"^ 

p n'^rj-f 1.3 X 10-^° 4.5 x 10"^° 5.2 x 10"^° 

a;^7r°7r°7 1.7x10"^ (2.9-3.5) x 10"^ (3.5-4.6) xlO-^ 
cj n^rj-f 3.4 x lO"'^ 9.9 x 10"^ 9.7 x 10"^ 

Table 2: Scalar (LcrM) and vector meson (VMD) exchange contributions to 
V P^P^^ decays. 

(j) TT^ri'j and p — > 7r°7r°7 can be used to extract relevant information on 
the properties of the /o(980), ao(980) and a (600) scalar mesons, respectively, 
while cu — > 7r°7r°7 and (p, a;) n^rj'j cannot be used due to the smallness 
of the scalar contribution, (j) K^K^^i could also be used in the near fu- 
ture for extracting the properties of the /o and resonances. The following 
comments are for each process in particular. For the 4> ~^ tt^tt'^j decay, we 
have seen that the /o exchange dominates whereas the contributions of the a 
and intermediate vectors are of 10% each. The dependence of our prediction 
on the /o mass and the scalar mixing angle could be used to obtain relevant 
information on these two parameters. We have also provided an explana- 
tion for the suppression of the a contribution. For — > n^rj'j, we have a 
parameter-free prediction which is dominated by the gq exchange and where 
the vector exchange is only of 5%. For the not yet measured — > K^K^^ 
decay, wc have shown that the /o contributes more strongly than the oq while 
other contributions are negligible. Our prediction seems to be very sensitive 
to the scalar mixing angle. A measurement of this process at DAFN&2 
would be welcome and can serve as an additional test of the whole approach. 
It is also confirmed that this process does not pose a background problem 
for testing CP- violation at DA$NE. Concerning the p — > 7r°7r°7 decay, our 
analysis predicts an important contribution from the a resonance which ac- 
counts for around half of the total signal. A comparison with experimental 
data seems to prefer a low mass and moderately narrow cr (600). Further 
experimental improvement on this decay is required before drawing definite 
conclusions. On the contrary, the scalar contribution to the uo 7r°7r°7 de- 
cay is negligible and the process is dominated by the intermediate p meson 
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exchange. However, this vector contribution complemented with p-uj mixing 
effects and the use of a momentum dependent p width is not enough to reach 
the large experimental branching ratio which remains unexplained. Finally, 
the decays (p, u) n^ri'y are not very appealing since the predicted size 
of their branching ratios is far from being observed in the near future. In 
conclusion, higher accuracy data and more refined theoretical analyses would 
contribute decisively to clarify the sector of the lowest lying scalar states. In 
this line, our work aims to be one step forward. 
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A Appendix 

A.l Scalar amplitudes 

K+K- K°K° 

The K+K- K^K^ amplitude in the LaM is 

jyhaM _ gaK+K-gaKORo 9 f()K+ R- 9 foKORO 

■^K+K-^KORO - 9K+K-KOR0 o - 771^ S - rn^ 

" ft) 

9aoK+K-9aoKORO ^a^KTKO(O) 



(31) 



where the SPP couplings in the isospin limit (see Refs. |33] for details) 
are 



9aK+K- = 9aK0R0 = 9aKR = -^^J^ [COS (ps - V2sm (ps) , 

2 2 

9hK+K- = 9foKORo = 9foKR = ^^^%^(sin05 + ^2 COS 0s) , (32) 
9aoK+K ~9aoKORO - ^ - 2f]^ , 

and (ps is the scalar mixing angle in the quark-flavour basis defined as 

a = cos (j)scrg - sin^scr^ , /o = sin^scXg + cos^^o-^ , (33) 
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with cTq = -^{uu + dd) and cr^ = ss. The gK+K-RORo couphng is fixed by 
imposing that the A^^.^j^of^o vanishes in the soft-kaon hmit (either p — >• 
or p' 0) implying s = t = u = m\. The resulting expression, 

2 2 (7^ - 

has been explicitly checked within the LaM. In consequence, 
.L<xM _ 91kr s-ml g%KR s-m\ 

'^R+ R-^R^R'^ 2 9 9~'~2 2 2 

— mj. s — — m^^^ s — mj^ 

2 

, gaoR+R-gaoRORO S - m]^ ^ ^a±/fT/fO(o) t - m\ 

+ 2 9 9~'~2 9J.9' y"^^} 

m\r — TTii s — mi ml- — mj, t — mi 
and utilizing the definite forms of the couplings in Eq. ( l32l) . one obtains 



2 

^JR 



^^^^^(c0,-v/2s0,)^ 
s — mi 



mi - m^„ ^ ^ X 9 - m^ 

__f^ ^ S05 + ' - — f 

s-mj^ s- ml^ 

t-m\m\- ml^ 



VI t - ml 



(36) 



In order to implement this amplitude into -^^^^o/fo^ preserving the chiral- 
loop prediction for infinite scalar masses, the former has to become first 
only s-dependent. As explained in Sect. El this is achieved replacing the t- 
dependent part of Eq. (136!) by the result of subtracting from the chiral-loop 
amplitude fill I) the infinite mass limit of the s-channel scalar contributions. 
In doing so, one finally gets the scalar amplitude 0141) where the naive scalar 
poles have been already substituted by more appropriate scalar propagators, 
namely a Breit-Wigner for the a and complete one-loop propagators for the 
/o and Oo (see Sect. I A. 2]) . The amplitude f|T^ . now only s-dependent, is then 



incorporated in -Ah^H^ofyo to give the result shown in Eq. ffT^ . 



^(I>->RORo^ 
-0^0 

The K^K^ vr'^Tr^ amplitude in the LaM is 



K+K- TT-TT- 



/(LcrM _ gaR+R-g<T-KOn° g foR+ R- g fo-K^-K^ 
^X+X-^ttOttO — gR+R-TvOwO 5 2 

s — mi s — mi 

JO 

~9l^R±-K"(': 2 ~^ 2) ' ^'^^^ 

\t — mi u — mij 
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where the required SPP couphngs, besides those of Eq. flH^ . are 



COSTS' 



9 9 

int-—m, 
^ to 



sm (ps 
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Notice that g^^K^-j^o can be written in two different equivalent forms. Using 
the soft-pion hmit (p^ — > 0), which imphes s = and t = u 



m\, the 



couphng gK+K-wOnO is fixed to 

gaK+K-gcTTTOnO g fQK+ R- g foirOnO 



gK+K-Tz^TrO 



ml — ml 



+ 



ml — m^^ 



2 9 ' 



(39) 



and the amphtude to 



A 



Lo-M 

K+K--*-kOttO 



gaK+K-ganOTrO S - ml ^ g foK+ R- g fonOnO S - ml 



ml — ml s 



mi 



ml — m'j^ 



s — m 



fo 



9l^K±7rO ft-ml ^u-m 



mj^ — \ t — m 



u — mt 



(40) 



that after utihzing the exphcit form of the couphngs in Eqs. 
becomes 



and 



/I Lo-M 



S ~ mt, 



K 



s — ml 



COS</>5(^COS (ps 



V2 



sm (ps) 



m^-mj^ . 



m 



sm (i)s[sm(f)s 



+ V2 



cos (ps) 



fo 



K 



(t - m\) 



m 



K 



m 



t — mi 



(41) 



One gets the s-dependent scalar amplitude in Eq. (1161) replacing the t- and 
u-dependent parts of the expression above by the corresponding chiral-loop 
amplitude ([H]) minus the infinite mass limit of the s-channel scalar contribu- 
tions in Eq. ( |4T1) . 



K+K- 

The K^K^ 71^1] amplitude in the LcjM is 



/I Lo-M 



gK+K-TvOri 



gaoK+K-gao 



s — m 



+ 



1 



t — mi u — mi 



(42) 
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where the required SPP couphngs are 

2 2 2 2 

— m — j- 

gaonOr, = 7 COS 0p , g^^K±l^ = TTJ. (cOS 0p - V 2 sill 0p) , (43) 

and (pp is the rj-r]' mixing angle in the quark-flavour basis defined as 

fj = cos (j)pr]q - sin (f)pT]s , t]' = sin (f)pr]q + cos (f)pT]s , (44) 

with rig = -^{uu + dd) and rjs = ss. Using the soft-pion limit (p^ — > 0), 

which implies in this case s = and t = u = mj^, the coupling QK+K-TrOri 
is fixed to 

9K+K-TT0r) — 2 2 ^^2 ' V^"-*/ 



and the amplitude to 



A 



2 

LaM ^ 9aoK+K-9ao-K°T} S - TU^ 

K^K~ — ^■Tr*^?? 9 9 9 



Using the explicit expressions for the couplings in Eqs. fl32|) . fl38l) and (H3l) . 
the former amplitude is written as 



+ 



fflO 

(t - m^)-^^ ^ + (« - m^)^^ ^ 

r — mf u — mi 



x-(cos0p — -\/2sin0p) I 



(47) 



Again, one gets the s-dependent scalar amplitude in Eq. ( fTTj) replacing the t- 
and M-dependent parts of the expression above by the corresponding chiral- 
loop amplitude ([9]) minus the infinite mass limit of the s-channel scalar con- 
tributions in Eq. fHTj) . 



_+_— , _o_o 



The TT^TT ^ TT TT amplitude in the LcxM is 

xLcrM „ 5'cr7r+7r-5'o-7r07rO 9 foir+TT- 9 foTT^TT^ /.o\ 
"A + 0^0 — 9n+n~n0n0 5 o ' 1^^°/ 

s — s — mi 

<^ JO 
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where the required SPP couphngs are 

2 2 

gan+n- = Qu-kO-kO = "/^ " COS 0g 



(49) 

9f07T+7T- = gfoTrOnO = Sm 05 , 



and the soft-pion hmit (p^r 0), which imphes s = t = u = m^, fixes the 
couphng g^+^-^o^^o to 

gn+TT-nOnO = 5 5 ^ 5 o " ^'^^) 

mi — mi mi — mi 

IT a T JO 

Hence, the amphtude is written as 

.LaU _ 9a7T+7T-9cTn0n0 S - ml ^'/ott+tt- ^'/pTrOTrO S - ml 

m^ — m^ s — m^ m^ — mj^ s — mj^ 

which turns out to be 

Al%-.^„„, = ( + =in^ I , (52) 

once the exphcit expressions in Eq. ( H9i) are used. As seen, this amphtude 
is already s-dependent and coincides with its chiral-loop counterpart ( fTOl) in 
the infinite scalar mass limit. Hence, it can be directly plugged into A^'^q^o^ 
to get the total scalar contribution to this process. 



A. 2 Complete one-loop propagators 

The complete one-loop propagators for the /o and scalar resonances are 
defined as 

D{s) = s-m\ + ReH(s) - ReH(m^) + ilmH(s) , (53) 

where m^ is the renormalized mass of the scalar meson and H(s) is the 
one-particle irreducible two-point function. ReH(m|.) is introduced to regu- 
larize the divergent behaviour of H(s). The propagator so defined is well be- 
haved when a threshold is approached from below, thus improving the usual 
Breit-Wigner prescription not particularly suited for spinless resonances (see 
Ref. [IB] for details). However, for the a meson we choose a simple Breit- 
Wigner propagator since the values of mass and width used in our analysis 
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are obtained from experimental data applying this parametrization. The real 
and imaginary parts of n(s) for the /o ar^ {Ris) = Ren(s), I{s) = Imn(s)) 



^fpKK 
167r2 



2 - log (i±|^ j - 2^^ arctan j 6 
2 - (3k log (i^lf ) 6;, - 2^K arctan (^) 



(54) 



I), 



where = a/1 — 4m^ / s for i = (vr, K), (3i = \/^mll s — 1, 9j = 9(s — 4m^^ , 
and G)j = 0(4m| — s). The couplings of the /o to pions and kaons are written 
in the isospin limit so g\^^ = lgj^^+^- and ^J^^^ = 2^J^^+^_. For the oq 
the real and imaginary parts of the two-point function are 



R(^) = '-h^[^-PK log (i^) 6;, - 2Pk arctan (^) 6;, 



167r2 



2-^logfe -/?^X-log 



-2^+^/?- arctan (||) 6., + ^+4" log (||±|| ) 6., 



(55) 



where = a/1 - (m^ ± rriny/s, (3^^ = ^J {m^ ± rrir^f j s - 1, 0^^ = e[s 
(m^ + m^)^], 6^^ = 6[s - (m^r - m^)^] x 6[(m7r + m.^)^ - s], and 6,^77 
e[(m^ - m^)2 - s]. 



A. 3 Invariant mass distribution 

The final P^P^ invariant mass distribution for the V P'^P^^ processes i^ 

dViy ^ P°P°7) ^ ^Fl^m ^ c^rvMD ^ ciFint ^gg^ 
dnipopoi dmpopoi dmpopoi dmpopoi ' 

® In our analysis, we work in the isospin limit and therefore the mass difference between 
and is not taken into account for the KK threshold. 

^ In terms of the photon energy, — {niyo — mpopo/)/(2myo), the photonic spectrum 
is written as dT/dEy = (myo/mpopoi) x dV / dmpapor . 



34 



where the LaM term (the scalar or scalar contribution) i a^l 



(57) 

X CLmI^(^?.opo,)I ^ P°P°)l.mP , 

with (|, l) for the case or not of identical pseudoscalar states in the final 
state, [g, Qs) for the case or not of the meson as the initial decaying vector, 
rrip the mass of the charged pseudoscalar inside the loop (m^+ for p — >• n^n^'j, 
mK+ otherwise), Clo-m the coefficients defined as 



Ct7M'" = -^20--"^^' = -V2C[-:^°^ = 1 , (58) 

and A{P~^P~ — > P^P'^)LaM the four-pseudoscalar amplitudes taken from 
Eq. (HM and Eqs. ( fTGUTHl) . The VMD term (the vector or background con- 
tribution) and the contribution resulting from the interference of both am- 
plitudes are given by 



'^r[vMD,int] _ /I 1 mpOpO, / ^ 

dmpOpOl V2' / 2567r3 m^o \ 



"pOpO/ 
— 

™ n 



1 

rfx A [vMD.int] ("^popo' , a;) 



1 



For the generic process V — > p'^p'^j the notation is fixed to define mpo as the mass 
of the first pseudoscalar in the reaction and mpo/ the mass of the second one. 
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where one explicitly has 

^VMD("^popo/, a^) = I J2po\ l«^VMDp = I ^C'vMD^^ 

X (I [mpoTJipo, - 2myompompo, + (m^o - my fm\,o + m\^m^po, 
+2mv{rnyo +rhv - my,)m],om],o, - 2my{my + my,)w?po, 
-2m^(m^o - m^)(m^o - ml - m*^, 
+my{myo + my - rhy,y] x 



pO 



rriymyi 



\Dv{rh 



and 



[mpo^po/ - rhy,{myo - rhyYm^po - rhy{myo - rhy,Ym\,o, 

+2mymy,{mya - rhy - my,)(mpo + mp,,,) + myrhy, 
-myfhy^fhy + fhy,){rnyo - rhy - rhy,)] 

^uit{m?pQ pQ, , x) = fR-G ^p^j ^l<tM>^vMD 



_ I fr< e(g,gs) 

Lo-M 



X- 



rh'^irny+niy, 



+ (mpo,m 



with (P2 = m^, P'2 = rh*^,) 



m 



V 



m 



*2 

v 



2 I "^yO "•'pOpO/ 

lll'pO I o 



1 + 



'po- 



''•pO/ 



'pO pO/ 



+ 



"^yO "''pOpO/ 



pO/ I 2 
pO ~l~ lll'pO, 



-y _ "^pO-'"pO/ 



pOpO' 



1 + 



+ 1 + 



"■pOpOl 



''pO/ 



'pO pO/ 



mid + m PO, + 



pOpO/ 



^2 



m 



V ' 



PpoX 
PpoX 



and 



4m^o 



^pOpO/ ( 1 _|_ '"pO ™pO/ 
V '"pOpO, 



2 • 
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The invariant mass distribution for the processes with mpo = m^po, can 
be simphfied to 



4opoA' 1 1 _ 

dnipopo 2 1927rB 47r m*^ m^o V"^ / V '"pOpO (64) 

X Cl,^\L{mlopo)\'\AiP^P- ^ ^°^°)LaM| 

for the scalar contribution, and 



<^r[vMD,int] _ 1 1 "tpOpO f-^_ ™pOpO^^ A _ '^™pO 

drripopo ~ 2 2567r3 m^,o \ m'^Q J V "^pOp( 



where 



X (I [m|,o - 2m^o"^po + (rnyo + 4m^ - 2m^m^2,)m^o 
-2fhvrriyo{rnyo - 2fhy - my,)m],o - Amy{my + my,)mpo 

+4mymy,{myo - rhy - rhy,)m'\,Q + rrtymy, 

-mymy,{my + my,)(m^o - ^Ti^ - rfi*^>)] 



and 



r^^(m^+m^V2m^o)^-m^Q(mf,-m^o)^ 



13^ (m2,) 

+ (m^,y)^(mt?„r)]} , 



(65) 



(66) 



^int (^po po ; 3;) — g Re X^poi »4.L(tM^VMD 

x2Re {L(m^opo)^(P+P- ^ P^'P^aU (67) 
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with 



o2 ^2 



4m' 
pO pO 



= ml. + ^^v^ll^ ( 1 +^ /l - ^ ] , (68) 



'pOpO 



2m PO + ml 



''po "T "tyo mpopQ m 
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